USAAVLABS  TECHNICAL  REPORT  68-67 


C7J 

o 

'  > 

\  «y 

CO 

CO 

^rf 


u 


LIFTING  SURFACE  THEORY  AHD  TAIL  DOWNWASH 
CALCULATIONS  FOR  V/STOL  AIRCRAFT 
IN  TRANSITION  AND  CRUISE 


ir 

E.  S.  Liviiskj 

I.  I.  Tkiiiti 
P.  M.  Yifir 
C.  N.  Nillail 


Octibir  1961 

U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTIS,  VIRGINIA 


CONTRACT  DAAJ02-67-C-0059 

AIR  VEHICLE  CORPORATION 
SAN  OIEGO,  CALIFORNIA 


This  document  has  been  approved 
for  public  release  and  sale ;  its 
distribution  is  uni imited. 


x 


JAN  21  SB 


Disclaimers 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  authorized 
documents. 

When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the 
said  drawings,  specifications,  or  other  data  is  not  to  be  regarded  by 
implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission,  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 


Disposition  Instructions 


Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to  the 
originator. 


DEPARTMENT  OF  THE  ARMY 

U.  a.  ARMY  AVIATION  MATtMKL  LABORATORY* 
PONT  tuarw.  VMKMMA  IH04 


Thla  raport  was  prapared  by  the  Air  Vehicle  Corporation  under  the 
terms  of  Contract  DAAJ02-67-C-0059.  It  conaiata  of  the  development 
of  a  theoretical  method  and  resulting  computer  program  for 
calculating  the  wake  flow  field  characterlatlca  of  wing  and  propeller 
comblnatlona  during  tranaltlon  of  V/STOL  vehicles. 

Reaaonable  agreement  haa  been  obtained  in  comparison  of  the  theory 
and  experimental  data  at  low  to  moderate  wing  angles  of  attack. 
However,  insufficient  experimental  data  are  available  to  test  the 
validity  of  the  theory  at  high  wing  angles  of  attack. 

The  report  is  published  for  the  exchange  of  information  and  the 
stimulation  of  ldeaa. 


Project  1F125901A142 
Contract  DAAJ02-67-C-0059 
USAAVLABS  Technical  Report  68-67 
October  1968 


LIFTING  SURFACE  THEORY  AND  TAIL  DOWNWASH 
CALCULATIONS  FOR  V/STOL  AIRCRAFT 
IN  TRANSITION  AND  CRUISE 


Final  Report 
Report  No.  356 


by 


E.  S.  Levinsky,  H.  U.  Thommen, 
P.  M.  Yager,  and  C.  H.  Holland 


Prepared  by 

Air  Vehicle  Corporation 
San  Diego,  California 


for 

U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTIS,  VIRGINIA 


Thla  document  has  been  approved 
for  public  release  and  sale;  its 
distribution  is  unlimited. 


V 


1 


A  large-tilt-angle  lifting- surface  ry  is  developed  for  tilt- wing  and 
tilt-propeller  (or  rotor)  type  V/STOL.  aircraft.  The  method  is  based  up¬ 
on  an  inclined  actuator  disc  analysis  in  which  closed-form  solutions  are 
obtained  for  the  velocity  potential  inside  and  outside  a  single  fully  con¬ 
tracted  circular  slipstream  at  large  distances  behind  the  actuator  surface. 
Both  the  normal  velocity  and  the  nonlinear  pressure  boundary  conditions 
are  satisfied  exactly  across  the  slipstream  interface  for  a  single  slip¬ 
stream  at  arbitrary  inclination  angles  to  the  free  stream. 

The  inclined  actuator  disc  analysis  is  combined  with  a  discrete-vortex 
Weissinger-type  lifting  surface  theory  for  application  to  wing-propeller 
combinations  at  arbitrary  wing  angle  of  attack,  propeller  tilt  angle,  and 
thrust  coefficient.  Wing-induced  modifications  to  the  boundary  conditions 
across  the  slipstream  interface  are  made  according  to  the  approach  of 
Ribner  and  Ellis.  The  theory  is  generalized  for  configurations  with  one, 
two,  or  four  slipstreams,  and  effects  of  slipstream  rotation  are  intro¬ 
duced  in  all  but  the  single- slipstream  cases. 

Sample  calculations  showing  a  comparison  with  previous  test  data  are 
presented.  Agreement  between  theory  and  experiment  is  shown  to  be 
satisfactory  for  small  slipstream  inclination  angles.  However,  at  large 
tilt  angles  the  theory  (with  an  undeformed,  but  displaced,  slipstream  and 
wake)  tends  to  predict  significantly  lower  downwash  angles  in  the  tail 
region  than  observed  from  a  single  set  of  test  data,  possibly  due  to  slip¬ 
stream  deformation  and  wake  roll-up  which  reduce  the  downward  wake 
displacement.  Use  of  only  one-half  the  calculated  wake  displacement 
gave  improved  agreement  with  downwash  angle  data  at  these  conditions. 
However,  insufficient  downwash  angle  data  are  available  for  making  a 
general  evaluation  of  the  theory  at  large  slipstream  angles. 

Extensive  digital  computer  results  are  given  in  chart  form,  showing  span 
loading,  downwash  angle  c,  stability  parameter  d«/da,  and  dynamic 
pressure  inside  as  well  as  outside  the  propeller  (rotor)  slipstream  at  ar¬ 
bitrary  points  behind  the  wing.  These  charts  are  presented  for  v/STOL 
configurations  with  two  and  four  slipstreams  and  at  flight  conditions  rang¬ 
ing  from  hover  to  cruise. 
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INTRODUCTION 


Considerable  research  effort  has  been  expended  in  the  past  pertaining  to 
the  interaction  between  a  wing  and  slipstream.  Most  noteworthy  is  the 
pioneering  contribution  of  Koning  (reference  1),  who  formulated  the  proper 
boundary  conditions  across  the  interface  of  an  uninclined  slipstream. 
Koning  matched  both  pressure  and  flow  angle  on  each  side  of  the  boundary 
within  the  limitations  of  linearized  theory  and  used  the  Prandtl  lifting  line 
theory  to  predict  slipstream  interference  effects  at  cruise  speeds.  Kon¬ 
ing1  s  treatment  was  extended  to  a  wider  range  of  forward  speeds  by 
Glauert  (reference  2)  and  by  Franke  and  Weinig  (reference  3).  The  latter 
authors  also  improved  and  generalized  the  theory  to  include  effects  of 
small  propeller  inclination  and  slipstream  rotation.  Nevertheless,  com¬ 
parison  with  test  data  indicated  that  the  theory  overpredicted  the  additional 
wing  lift  inside  the  slipstream  by  a  substantial  margin  (reference  4). 

Partly  because  of  this  failure  of  the  lifting  line  approach,  which  was  as¬ 
sumed  due  to  the  low  aspect  ratio  of  the  wing  segment  inside  the  slip¬ 
stream,  Graham,  et  al.  (reference  5),  supplemented  the  lifting  line  theory 
with  lifting  surface  (Weissinger)  and  slender  body  theories,  and  found 
improved  agreement  with  test  data.  The  Weissinger  lifting  surface  the¬ 
ory  was  subsequently  generalized  in  a  series  of  papers  by  Ribner  and 
Ellis  (references  6  through  8)  to  apply  to  multiple  uninclined  slipstreams 
of  arbitrary  cross  section.  The  hitter  authors  introduced  horseshoe  vor¬ 
tex  elements  of  unknown  strength  along  the  wing  l/4-chord  line,  and  re¬ 
quired  that  the  wing  boundary  condition  be  satisfied  at  the  3/ 4-chord  line 
in  accordance  with  the  standard  Weissinger  approach.  By  considering  a 
reduced  velocity  potential  inside  the  slipstream,  they  were  able  to  satisfy 
both  the  pressure  and  flow  angle  boundary  conditions  across  the  slip¬ 
stream^)  through  introduction  of  an  additional  set  of  unknown  horseshoe 
vortex  elements  around  the  slipstream  periphery.  Calculations  carried 
out  for  a  wing  and  a  single  uninclined  slipstream  showed  good  agreement 
with  the  test  data  of  Brenckmann  (reference  9). 

It  is  clear  that  the  theory  of  wing-propeller  interaction  must  be  extended 
to  include  effects  of  large  slipstream  inclination,  if  the  theory  is  to  be 
applicable  for  tilt-wing  or  tilt-rotor  type  V/STOL  aircraft.  Previous 
attempts  to  treat  inclined  slipstreams  have  been  unsuccessful,  partly  be¬ 
cause  of  the  use  of  a  solid  cylinder  approximation  for  the  inclined  slip¬ 
stream,  e.g.  ,  references  5  and  8.  Although  the  solid  cylinder  produces  a 
disturbed  flow  field  outside  the  slipstream  which  satisfies  the  normal 
velocity  boundary  condition,  the  nonlinear  pressure  boundary  condition  is 
not  satisfied.  The  rectification  of  this  inadequacy  and  the  development  of 
a  wing- slipstream  interaction  theory  applicable  to  V/STOL  technology 
are  the  motivations  for  the  present  study. 

In  the  present  report  a  new  inclined  actuator  disc  theory,  which  satisfies 
both  the  normal  velocity  and  nonlinear  pressure  boundary  conditions 
across  the  slipstream  interface,  is  first  derived.  The  inclined  actuator 
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disc  theory  is  subsequently  combined  with  the  Ribner-Ellis  li/ting  surface 
method  to  treat  tilt-rotor  or  propeller-wing  combinations  at  large  tilt 
angles  and  at  forward  speeds  from  hover  to  cruise.  Effects  of  slipstream 
rotation  are  included  by  adopting  a  viscous  core  vortex  model  inside  the 
slipstream  for  the  swirl  velocities.  The  major  assumptions  in  the  re¬ 
sulting  generalized  analysis  are  that  the  slipstream  is  fully  contracted,  is 
of  basically  circular  cross  section,  and  is  of  constant  total  head  (but  with 
nonuniform  velocities  and  static  pressures)  in  the  region  of  the  wing. 
Unseparated  flow  is  assumed  over  the  wing  panels,  both  inside  and  out¬ 
side  the  slipstream,  and  effects  of  wake  roll-up  and  slipstream  distortion 
are  assumed  to  be  negligible. 

The  theory  is  compared  wherever  possible  with  available  experimental 
data.  In  addition,  an  extensive  series  of  calculations  is  carried  out  to 
provide  theoretical  data  on  the  span  loading  and  on  the  downwash  angles 
and  dynamic  pressures  in  the  wake  and  slipstream  for  representative  two- 
and  four-slipstream  V/STOL  configurations  and  flight  conditions.  These 
results  are  presented  in  chart  form  for  the  rapid  estimation  of  the  stabil¬ 
ity  and  control  characteristics  and  of  wing  structural  requirements.  It  is 
recognized  that  additional  experimen.al  verification  of  the  theory  is  re¬ 
quired,  especially  for  large  tilt  angle  conditions  for  which  detailed  wake 
and  slipstream  surveys  are  largely  unavailable. 


INC  LINED  ACTUATOR  DISC  THEORY 


BOUNDARY  CONDITIONS  AND  GENERAL  SOLUTION 

The  incompressible  flow  field  created  by  an  actuator  disc  with  its  axis  in- 


(Figure  1).  The  air  passing  through  the  actuator  disc  is  assumed  to 

experience  a  uniform  increase  in  total  and  static  pressure  of  magnitude 

AH  ,  so  that  the  total  pressure  inside  the  slipstream  has  the  constant 

value  H  =  H  +  AH,  where  H  is  the  total  pressure  outside  the  slip- 
s  o  o  r  r 

stream.  The  resultant  force  on  the  actuator  disc  is  assumed  to  act  nor¬ 
mal  to  the  disc  surface,  and  is  of  magnitude  A^  AH.  In  this  respect 

the  actuator  disc  differs  from  an  inclined  propeller  or  rotor  which  devel¬ 
ops  forces  in  the  plane  of  rotation.  This  in-plane  force  is  usually  small 
compared  to  the  normal  force,  and  is  neglected  in  the  present  formulation. 

We  next  assume  that  this  flow  field  can  be  approximated  far  behind  the 
actuator  disc  by  a  distribution  of  suitable  singularities  around  the  surface 
of  a  semi-infinite  circular  tube  (Figure  1).  The  axis  of  the  tube  is  inclined 
to  the  free  stream  by  the  angle  6.  This  is  the  angle  at  which  the  slip¬ 
stream  would  be  inclined  far  behind  the  disc  if  not  for  distortion.  The 
radius  of  the  tube  is  taken  as  that  of  the  fully  contracted  slipstream  r  =  R, 
and  the  boundary  conditions  which  must  hold  across  a  fluid  interface  are 
satisfied  on  the  surface  of  the  tube. 


A  slipstream  coordinate  system  x,  y,  z  (or  equivalently  x,  r,  0)  with  the 
x-axis  coincident  with  the  tube  axis  is  introduced  for  convenience.  De¬ 
noting  by  V  ,  V ^  the  velocities  outside  and  inside  the  slipstream,  re¬ 
spectively,  the  boundary  condition  on  the  velocity  components  normal  to 
the  slipstream  boundary  is 


e  • 
x 


V 

o 


for  r  =  R 


(1) 


where  e^,  er>  e^  are  unit  vectors  in  the  designated  directions  (Figure  1). 

Equation  (1)  signifies  that  the  flow  angles  normal  to  the  cylindrical  sur¬ 
face  are  the  same  across  the  interface,  but  allows  a  discontinuity  in  the 
velocity  components  in  the  piine  of  the  surface. 

The  axial  velocity  components  of  Vq  and  Vg  are  designated  by  uq 

and  ex  u^,  respectively,  and  may  be  regarded  as  constant  at  large 

distances  from  the  actuator  disc.  However,  the  cross-flow  velocity 
components  are  nonuniform,  even  for  large  x,  and  lead  to  deformation 
of  the  originally  circular  slipstream. 
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A  second  boundary  condition  is  that  the  static  pressures  pQ(r,  0)  and 

PB(r,  0)  outside  and  inside  the  slipstream,  respectively,  are  identical 

on  the  tube  surface.  It  will  prove  convenient  in  the  analysis  to  change  the 
flow  into  one  of  constant  total  head  by  subtracting  AH  from  the  pressure 
ps  inside  the  slipstream.  The  pressure  boundary  condition  is  then 

Po(R,  0)  -  pa(R,0)  =  AH  (2) 


The  above  flow  field  may  be  obtained  by  suitably  distributing  singularities 
on  the  surface  of  the  cylindrical  tube.  First,  consider  a  distribution  of 

vortices  y  =  eQ  ^  j  pointed  in  the  azimuthal  direction  with  strength  T 

independent  of  0(ring  vortices).  Restricting  ourselves  to  large  distances 
behind  the  actuator  disc,  the  ring  vortex  distribution  may  be  regarded  as 
infinitely  long  and  of  uniform  strength.  It  will  therefore  induce  only  axial 
velocity  components  inside  the  tube,  and  will  have  no  influence  outside 
the  tube.  The  derivative  y  is  chosen  so  that  the  induced  velocity  V 
-.  V 

inside  the  tube  is  e  (u  -  u  )  ,  viz.  , 

x  s  o 


lim  V  =  e  y 
y  x  ' 

X-*-CO  ' 


e  ( u  - 
x  s 


%> 


lim  V  =  0 
x  -»  oc  Y 


r  <  R 

(3) 

r  >  R 


Next,  a  distribution  of  ”ortices  and  sources  and  sinks  with  axes  parallel 
to  x  is  introduced.  These  singularities  will  vary  in  strength  as  a  func¬ 
tion  of  0,  but  can  be  assumed  independent  of  x  when  x-*-®.  In  this 
limit,  and  using  symmetry  conditions  about  0  =  t  tt/2  ,  the  most  general 
form  of  the  resulting  two-dimensional  velocity  perturbation  potential  may 
be  written  as 

£  RCn(r/R)n 

=  <P  g  =  "  /  n'  sinn0  ,  r  <  R  (4) 

n=l ,  3,  5 


V  RDn(r/R)"n 

<p  =  <p  =  )  - -  sin  n  G  , 


n= 1 .  3, 5 

Thus,  the  resultant  velocities  for  x-*®  become 


V  =  V  +  e  (u  -u)+V<p 
s  «  x  s  o  s 


e  u  +  e  w  +  V  <p 
x  s  z  «  s 


r  >  R  (5) 


r  <  R  (6) 


V 

o 


+  V?0 


w  +  v  ffl 
oo  v  r  o 


r  >  R  (7) 
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where  the  undisturbed  free-stream  velocity 
components 


V  has  been  resolved  into 
00 


V  =  e  u  4  e  w 
ao  x  o  z  oo 


with 


u  =  V  cos  6  and  w  =  V  sin  6 

O  00  00  CO 


EVALUATION  OF  THE  UNKNOWN  COEFFICIENTS 


The  coefficients  C  ,  D  ,  which  must  be  known  to  obtain  the  velocity 
n  n  7 

field,  are  determined  from  the  boundary  conditions  on  the  slipstream 

interface. 

Making  use  of  Eqs.  (6)  and  (7),  the  normal  velocity  condition  [  Eq.  (1)] 
becomes 


8^o/3r  -  p.8^  /3r  *  -  (1  -  w^  sin8,  r  =  R  (8) 


where  u.  =  u  /  u  . 

^  o'  s 


Introducing  Eqs.  (4)  and  (j),  and  satisfying  Eq.  (  8)  for  each  term  in 
sin  n9  ,  gives  the  following  conditions  on  the  coefficients: 


Dj  *  kCj  =  (1  -  \)  wB 

D  -  jiC  , 
n  n 


(9) 

n  '  i 


Utilizing  the  steady  form  of  the  incompressible  Bernoulli  equation,  the 
pressure  condition  [  Eq.  (2)]  may  be  similarly  expressed  as 


2 


AH 

P 


2  r- 

-  u  4  e 


w  4  V 

z  e 


'J-f 


e  w  4  V  w 

7  00  T  O 


r  =  R 


(10) 


where  the  density  p  is  faken  as  the  same  inside  and  outside  the  slip¬ 
stream.  Since  7®  and  v©  are  both  functions  of  8,  whereas  the  re- 

s  o 

maining  quantities  in  Eq.  (10)  are  independent  of  8,  it  is  required  that 


or,  equivalently, 


2  AH  2 
=  u 

P  8 


AH 

9co 


1  -  u  2  c 

— cos  o  > 

V- 


r  =  R 


r  =  R  (11) 
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in  order  that  Eq.  (10)  be  satisfied.  Using  complex  n^'  *tion,  the  pressure 
condition  becomes 

liw  +  V  I  =  I  i  w  +  Vfli  I  ,  r=R 

1  CD  s  I  i  oo  *o  1 

which  may  be  written,  without  loss  of  generality,  as 

e*  '  ^  (w  -  i  V(?  )  =  w  -  i  V  <p  ,  r  =  R  (12) 

'  CO  S  *0  r  o 

where  £  is  an  unknown  function  of  9  to  be  determined,  and 
V  <p  -  9^)/0y  +  i  (d<p/dz).  Substituting  Eqs.  (4),  (5),  and  (9)  into  Eq.  (12) 
for  r  =  R,  and  observing  that  all  coefficients  Cn,  Dn  are  real,  it  is^ 
readily  determined  that  £,  must  be  either  0  or  it.  Taking  first  £,  =  0  , 

:  o 

and  equating  coefficients  of  like  powers  of  e  on  both  sides  of  Eq.  (12), 
gives 


C3  =  -(1-H)wm  (13) 

C!=-pC,  n  2:  3 
n+2  r  n 


Complex  Perturbation  Potential 

In  terms  of  the  complex  potentials  F  (  =  <p  4  i  )  and  F  (=  <p  ■+  i 4*  )  in” 

r  r  s  s  s  o  o  o' 

side  and  outside  the  slipstream,  respectively,  Eqs.  (4)  and  (5)  yield 

dF  /d  Z 
s 

d  F  /d  Z 
o 

where  dF/dZ  =  (<p  r  cos  0  ■ 

a  v  i  9 

and  Z  =  r  e 

By  inserting  Eqs.  (9)  and  (13)  for  the  coefficients  C  ,  in  Eq.  (14),  it 
is  found  that  the  series  can  be  readily  summed.  The  resulting  closed - 

*  The  rejection  of  solutions  with  t  =  -rr  is  justified  in  the  next  section. 


2  C„<z/R> 
- 1 ,  2 

‘I 


n=l ,  3, . 


(14) 


D 


n=  1 ,  3,  . 


(Z/R) 


n-f  1 


<p  sin  0)  -  i  (<p  r  sin  0  4  ^  <p  cos  0) 
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form  expressions  may  be  integrated  with  respect  to  Z  ,  giving 


In  the  limit  p-»l,  which  corresponds  to  an  actuator  disc  of  vanishing 
strength,  AH  +  0  and  Fg  and  Fq  vanish,  indicating  that  the  flow  field 

is  undisturbed  in  this  limit.  Re-evaluation  of  the  coefficients  for  £  =  tt 
yields  expressions  similar  in  form  to  Eqs.  (15)  and  (16),  except  that  the 
factor  (l-p)is  replaced  by  (1  +  p).  These  solutions  must  be  discarded, 
because  they  fail  to  satisfy  the  physical  flow  in  the  limit  p-M. 

In  the  limit  p-*-0,  corresponding  either  to  an  inclined  actuator  disc  of  in¬ 
finite  strength  (AH*'®)  or  at  zero  forward  velocity,  expansion  of  the 
logarithmic  term  in  Eq.  (16)  shows  that  Fq  becomes  the  complex  poten¬ 
tial  produced  by  a  doublet  at  the  origin.  Upon  superposition  of  the  com¬ 
ponent  w^  of  the  free  stream  velocity  parallel  to  the  z-axis,  the  result¬ 
ant  outside  flow  velocity  becomes  the  same  as  if  the  slipstream  boundary 
were  a  solid  surface.  The  solid  cylinder  approximation  for  an  inclined 
slipstream,  which  was  discussed  in  the  Introduction,  is  found  to  hold  only 
in  this  limit. 


At  intermediate  values  of  p  ,  the  complex  perturbation  potential  F0  is 
recognized  as  that  due  to  a  source  of  strength  w  (1  -  p)  w^  R/ p^^  located 
at  Z  =  -iRp^/^  and  to  a  sink  of  equal  strength  located  at  Z  =  i  R  p*'  ^ . 
Similarly,  Fg  is  that  of  a  source  and  sink  at  Z  =  t  i  R/p  together 

with  a  parallel  flow  of  velocity  ( 1  -  p)  w^/  p  in  the  positive  Z  direction. 


Velocity  Field 


Having  obtained  the  complex  potential  in  the  Trefftz  plane  x  =  <*> ,  the 
overall  velocity  field  is  readily  evaluated  from  Eqs.  (6)  and  (7).  The 
velocity  potentials  <p  and  <p  q  are  found  by  extracting  the  real  parts  of 

Fg  and  Fq  ,  respectively.  Thus,  in  rectangular  coordinates, 


(17) 


1 
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(1  -  u)  „  ,  |  v^  +  +  2Rz 


(18) 


Taking  the  gradient  of  tne  above  expressions  and  substituting  into  Eqs.  (6) 
and  (7),  the  overall  velocity  components  u,  v,  w  in  the  x,  y,  z  direc¬ 
tions  at  any  arbitrary  point  in  the  flow  field  are  found  to  be 


u  =  V  cos  6 


2Al  -  Wy*. 


M.u-r,'/2.' 


R  V  sin  6 


w  =  J 
o  * 


1  + 


(1  -u)R2(y2-s2tnR2) 


Tv24 


2  2  2 
for  y  +  z  >  R  , 


(z  -  Rp 

and 


*/2,  iry2+u+R111/2,' 


V  sin  6 
00 


(19) 


u  =  V  cos  b/\i 

8  CO  '  n 


2,1/2 

py  +  (  p  z  -  R) 


py  Mp 


>/2,* 


R) 


w 


1  - 


r,2+<r1/2*-R)2 


p  y  +  ( p 


>  R  V  sin  6 
'  oo 


R) 


sin  6 


2  2  2 
for  y  +  z  <  R  . 


(20) 


Equations  (19)  and  (20)  give  the  velocities  in  the  flow  field  at  x  =  oo  in 
terms  of  the  parameters  p  and  6.  Before  applying  these  expressions  to 
evaluating  the  downwash  angle  distribution  inside  as  well  as  outside  the 
slipstream,  it  will  prove  convenient  to  express  p  and  6  in  terms  of 
parameters  of  more  physical  significance.  This  will  be  accomplished  in 
the  next  section  by  relating  p  and  6  to  the  tilt  angle  and  thrust 

coefficient  T".  It  will  also  be  convenient  to  relate  the  vortex  tube  radius 
c 

R  to  the  physical  radius  of  the  disc  R  . 
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Forces  and  Power  on  an  Inclined  Actuator  Disc 

In  order  to  evaluate  the  forces  acting  on  the  disc,  a  cylindrical  control 
tube  of  radius  is  placed  around  the  disc  and  slipstream  as  shown  in 

Figure  2.  The  value  of  is  taken  sufficiently  large  so  that  the  pres¬ 

sure  is  ambient  around  the  cylindrical  surface  S^.  The  velocity  is  as¬ 
sumed  to  be  equal  to  V  on  the  forward  face  S.  at  x  =  -  as,  and  to  be 
given  by  Eqs.  (19)  and  (20)  on  the  rear  face  at  x  =  ®.  The  continuity 
equation  is  written  in  integral  form  as 


V  (-  ®)  •  d  S  4 


M 

s,  s, 

2  s  2  o 


V  (®)  •  dS  4 


I  7(Rc 


)  •  dS  =  0 


where  dS  has  the  direction  of  the  inward  normal  to  the  surface,  and  the 
subscripts  s  and  o  designate  the  portions  of  inside  and  outside  the 

slipstream,  respectively.  Using  Eqs.  (6)  and  (7),  the  net  mass  flux  into 
the  cylindrical  surface  becomes,  for  a  unit  density, 

f  V  (R  )  •  dS  =  ttR2  (u  -  u  ) 

.)  c  so 


and  carries  a  net  momentum  flux 


V  (V  •  dS)  =  itR  (u  -  u  )  V 

so® 


The  momentum  equation  is  expressed  in  the  form 

T  =  -  p  y  V  (V  •  dS)  -  y  pds“  (22) 

VS2+  ^3  ^l4  ^24  ^3 

where  T  is  the  thrust  force  normal  to  the  actuator  disc.  Taking  the 
component  of  T  in  the  x  direction  and  making  use  of  Eqs.  (6),  (7),  ( 1 1 X 
and  (21),  together  with  the  incompressible  form  of  Bernoulli's  equation, 
after  some  rearrangement  we  obtain 


9 


2 


X  =  irR2 


Here  we 


where 


/  ..  \  _  1  ~  H —  pw  tt  R 

PU8  (U8  ’  V  P  “ 


-  p  w 

'  n 


■o  +  -2  +  l,Got7,]'R 


L  V- 


(23) 


have  defined 


I  = 


9  <P  c 

8  * 


i 

Z  o 


_  9  <p 

i  r  __i 
=  “TT  j  9  z 

1T  R  % 


dS 


d  S 


'2  s 


;°  '  1 

b2o 


A2  I 


(V*Q)  dS 


G  = 
s 


(V  i  )  d  S 


’2  s 


is  a  modified  velocity  potential  such  that 

0*o/9z  =  (i/wa)  9<f>0/dz 
a^o/8y  =  (1/ W^  )  9^0/9y 
8^  g/8z  =  (m-/  wm  )  (9  <P  s  /  9z^  '  ^ 1  '  ^ 
9  y>  /  9y  =  ^/wco  }  9’’s//0y 


(24) 
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The  integrals  given  by  Eq.  (24),  as  evaluated  in  Appendix  I  ,  are 


!o  =Is  =  -(1"^ 


G  =  G  =  In 

O  S  L  p 


1 ±JL) 


Substituting  these  results  into  Eq.  (23)  gives 
X  =  ^Rzp  'us(us-uo)  +  |(1  +  R2) 

(25) 

Carrying  out  a  similar  evaluation  for  the  component  of  T  in  the  z  direc¬ 
tion  (the  pressure  forces  do  not  contribute  to  this  component)  gives 


Z  -  tt  R  p  ( 1  -  p)  uq 


(26) 


It  is  convenient  to  express  the  forces  on  the  actuator  disc  in  terms  of 
horizontal  and  vertical  force  coefficients  C^"  and  C^"  ,  based  on  the 

P  P 

propeller  area  and  the  dynamic  pressure  in  the  slipstream  at  ambient 

static  pressure.  Thus  we  define 

C = 


,  X  cos  6  -  Z  sin  6 
'H  ”  A  (q  +  AH) 

p  n(D 


„  ,,  _  X  sin  6  +  Z  cos  6 
CL  =  A  (q  +  AH) 

p  p 

Substituting  for  Z  and  X  through  Eqs.  (25)  and  (26),  and  using  Eq.  (11) 
for  AH  ,  after  some  simplification  we  obtain 


/A 

P  v  P 


V  =  br)  - 


|2(  1  -  p)  -  (  1  -p)  l  1  +  p  +  g(p)]  sin26|  cos  6 


x  -(f-) 

D  \  D  ' 


J^(l  -  p2)  cos26  +  p2  j 
|2  ( 1  -  p)  ( 1  +  p2)  -  (1  -  p)  [  1  +  p  +  g  (p)  ]  sin2 6  j-sin  6 


(27) 


T/i  Z\  2.  2 

(1-^i  )  COS  6  +  ^ 


(28) 
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where  g  ((a)  = 
ted  slipstream 


(j^H 

(vortex  tube). 


and  A  =  ttR^  is  the  area  of  the  fully  contrac- 

s 


Since  T  is  assumed  to  act  normal  to  the  plane  of  the  disc,  the  tilt  angle 
is  given  by 


a  =  tan 
P 


‘(V/v) 

'  D  D  / 


P  P- 

Introducing  the  thrust  coefficient  T"  ,  defined  as 


(29) 


T/A 

•p  '  i  —  P 


q  +  A  H 


y * (c-:)2] 


1/2 


the  ratio  of  the  propeller  area  to  slipstream  area  A  /A  becomes 

p  s 


A  /A  = 
P  3 


VA,)2+  (Ch'dA,/Ap)2J1/2/Tc 


(30) 


where  the  quantities  C^’’  A  / Ag  and  C^"  A  / Ag  are  given  by 

P  P  P  9 

Eqs.  (27)  and  (28). 

Equations  (11)  and  (27)  through  (30)  may  be  used  to  express  the 
flow  field  velocities  in  terms  of  the  given  physical  parameters  T  " ,  , 

and  A  instead  of  u,  6,  and  A  .  Once  T ",  or  ,  and  A  are  known, 
p  n  s  c  p  p 

(i,  6,  and  Ag  may  be  obtained  either  by  iterating  on  the  value  of  6,  or 

graphically  from  Figure  3. 

The  power  P  added  to  the  slipstream  and  to  the  outside  flow  by  the  actu¬ 
ator  disc  is  readily  obtained  by  evaluating  the  net  flux  of  energy  out  of  the 
cylindrical  control  surface  in  Figure  2.  Thus,  writing 


P  =  -  J  (V  V)  (V  •  dS)  -  J  pV-dS 
S 1 +  S2+S  3  S 1 +S2+S  3 


and  performing  a  similar  analysis  to  that  done  previously  for  the  com¬ 
ponents  of  "T  ,  we  obtain 
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(31) 


P  =  u  A  AH 
s  s 

The  effect  of  propeller  tilt  angle  on  power,  for  a  fixed  forward  velocity 
and  thrust  coefficient,  is  obtained  by  writing  Eq.  (31)  as 


Pla  ) 

. 

TrToT 


rp  1  t 

2,  .  c 

1/2 

A 

A^y 

p 

COS  0  T  J  f  I 

c 

» 1 

A 

1  1  c 

s  /  \ 

1  -  T" 

—  (0) 

c 

L  p  J 

(32) 


where  6  and  A(  / Ap  are  obtained  in  terms  of  a  from  Eq  (11)  and 
Eqs,  (27)  through  (30).  This  expression  has  been  evaluated  (Figure  4), 
and  is  compared  with  wind  tunnel  measurements  of  power  versus  tilt  angle 
in  the  following  section. 


COMPARISON  WITH  PROPELLER  TEST  DATA 


Because  of  the  many  assumptions  inherent  in  the  inclined  actuator  disc 
theory,  it  is  of  importance  that  the  theory  be  evaluated  by  comparison 
with  test  data  before  being  applied  to  the  wing-slipstream  interference 
problem.  Comparisons  of  the  power  and  average  downwash  angle  across 
the  slipstream  may  be  made  without  consideration  of  slipstream  rotation 
or  swirl.  However,  the  swirl  effects  must  be  added  to  the  theory  before 
detailed  downwash  survey  data  may  be  compared. 

A  comparison  of  predicted  and  measured  effects  of  propeller  tilt  angle  on 
power  is  shown  in  Figure  4.  Both  the  theory,  Eq.  (32),  and  the  test  data, 
reference  10,  show  that  the  power  required  to  produce  a  given  value  of 
T^'  decreases  with  increasing  till  angle,  with  the  power  becoming  inde¬ 
pendent  of  tilt  angle  as  T ^ '  ■*  1. 0.  The  data  scatter  and  differences  be¬ 
tween  theory  and  test  data  in  Figure  4  could  be  due  to  errors  in  estimating 
the  ideal  efficiency  of  the  inclined  propellers. 

The  downwash  angle  c  far  behind  the  propeller  may  be  found  from  the 
expression 


c  (y, z) =  -  tan 


- 1  f  w  cos  6  -  u  sin  6 
j  w  sin  6  4  u  cos  6 


(33) 


and  is  readily  determined  from  the  theory  using  Eqs.  (19)  and  (20).  As 
illustrated  in  Figure  5,  Eq.  (33)  predicts  a  nonuniform  downwash  distri¬ 
bution  inside  the  slipstream  and  a  varying  upwash  outside  the  slipstream, 
except  for  the  limiting  cases  T"  =  0  and  1.0.  The  calculated  downwash 


1  3 


angle  at  the  center  of  the  slipstream  has  been  plotted,  for  convenience, 
in  Figure  6. 

The  nonuniform  value  of  c  makes  comparison  with  test  data,  obtained  by 
means  of  a  tailplane  spanning  the  slipstream,  somewhat  difficult  to  inter¬ 
pret.  Nevertheless,  as  illustrated  in  Figure  7,  such  "average"  test  data 

(reference  11)  show  a  variation  in  the  ratio  of  c/ a  with  T  "  similar  to 

'pc 

that  predicted  by  theory,  at  least  for  the  small  values  of  a  for  which 
the  tests  were  conducted.  ^ 

A  simple  vortex  model  of  the  slipstream  has  been  used  to  calculate  the 
swirl  velocities  caused  by  the  propeller  rotation  (Appendix  II ) .  The  analy¬ 
sis  is  similar  to  that  by  Franke  and  Weinig  (reference  3)  in  that  a  central 
vortex  is  introduced  along  the  slipstream  axis  and  a  helical  vortex  is 
placed  around  the  slipstream  interface.  However,  in  order  to  keep  the 
swirl  velocities  finite  near  the  axis,  the  central  vortex  is  "softened"  by 
introducing  viscous  core  effects  according  to  Lamb  (reference  12).  The 
resulting  expression  (from  Appendix  II)  for  the  circumferential  swirl 
velocity  v^(r)  inside  the  slipstream  is 


(34) 


where  J  =  V  /2NR  is  the  advance  ratio,  v,  is  the  effective  turbulent 
<*>'  p  t 

viscosity,  ft  is  the  effective  starting  distance  ahead  of  the  origin  for  the 
viscous  core,  and  the  sign  is  positive  for  a  counterclockwise  sense  of 
rotation  when  viewed  from  the  rear.  No  swirl  velocities  are  predicted 
outside  the  fully  contracted  slipstream. 


Equation  (34)  is  compared  with  measurements  by  Stiiper  (reference  4)  in 
Figure  8.  The  test  data  were  obtained  at  several  stations  in  the  wake  of  a 
0.  49-ft. -diameter  pusher  propeller  aligned  with  the  free  stream  and  at  a 
thrust  coefficient  T"  ■=*  0.  38.  With  the  origin  x  =  0  located  at  the  blade 


trailing  edge,  the  best  fit  to  the  data  was  accomplished  by  choosing 
=  0.093  ft^/sec  anti  ft  *  0.08  ft  (approximately  twice  the  maximum 


blade  chord).  The  propeller  advance  ratio  was  0.  472,  whereas  the  free- 
stream  velocity  was  approximately  100  fps.  The  circulation  F  was 

14.  2  ft*"/  sec  as  computed  from  Eq.  (51)  of  Appendix  II. 


It  is  not  clear  how  will  scale  for  different  values  of  these  parameters, 

although,  as  discussed  by  Hall  (reference  13),  vt  scaling  with  either  T  or 

(Dl/  ^  has  been  considered.  Because  it  would  be  premature  to  use  such 
scaling  laws  indiscriminately,  and  because  of  the  lack  of  additional  test 
data  on  the  slipstream  vortex,  it  was  decided  to  use  =  0.  10  ftv  sec  for 

all  subsequent  calculations  and  design  charts  employing  swirl. 
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A  comparison  of  Eq.  (34)  with  test  data  from  reference  14  is  shown  in 
Figure  9.  The  test  data  were  obtained  at  a  single  station  located  1.43  ft. 
behind  the  propeller  (3.  25  ft.  diameter).  The  propeller  was  mounted  on  a 
flat-faced  cylindrical  nacelle  which  extended  rearward  to  the  measuring 
station.  This  may  have  been  the  cause  of  the  thick  viscous  core  which  is 
apparent  from  the  data.  In  order  to  achieve  correspondence  between  the 
theory  and  data,  it  was  necessary  to  place  the  effective  turbulent  origin 
28.  6  ft.  ahead  of  the  measuring  station  (20  times  the  actual  distance) 


when  using  a  nominal  value  of 


0.  10  ft2/ 


sec. 


Although  the  correla 


tion  shown  in  Figure  9  is  reasonable,  data  scatter  (due  possibly  to  wind 
tunnel  wall  interference  at  propeller  incidence)  and  the  scarcity  of  data 
points  make  the  comparison  somewhat  questionable. 


Additional  experimental  verification  of  the  inclined  actuator  disc  theory  is 
clearly  needed,  both  with  and  without  effects  of  slipstream  rotation.  This 
is  especially  true  for  conditions  of  large  tilt  angle  and  thrust  coefficient 
representative  of  V/STOL  transition,  for  which  little  experimental  data 
are  available.  In  the  next  section  the  inclined  actuator  disc  theory  is  used 
in  the  development  of  a  wing  lifting  surface  theory  with  inclined  slip¬ 
streams. 
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WING  LIFTING  SURFACE  THEORY  WITH 
INCLINED  SLIPSTREAMS 


A  method  will  be  developed  for  calculating  the  span  loading  and  wake  flow 
field,  inside  and  outside  the  slipstream,  for  a  wing  with  inclined  slip¬ 
streams  representative  of  tilt-wing  or  tilt-propeller  (rotor)  V/ STOL  air¬ 
craft.  The  procedure  will  be  applied  to  aircraft  with  one,  two,  or  four 
slipstreams.  For  reasons  of  computational  simplicity,  the  plane  of  the 
wing  will  be  assumed  to  pass  through  the  center  of  the  slipstream(s)  and  to 
bisect  the  slipstream  into  upper  and  lower  halves.  Also,  for  reasons  of 
simplicity,  the  formulation  will  be  restricted  to  configurations  with  sym¬ 
metrical  span  loadings  about  the  aircraft  centerline.  This  eliminates 
consideration  of  slipstream  swirl  for  single  slipstream  aircraft,  and  re¬ 
quires  that  the  propellers  be  rotating  in  opposite  directions  on  each  wing 
panel  for  the  two-  and  four-slipstream  configurations. 

Basically,  the  method  to  be  developed  can  be  considered  a  combination  of 
the  inclined  actuator  disc  theory  developed  in  the  previous  sections  and  the 
Ribner  and  Ellis  potential  theory  for  a  wing  in  an  uninclined  slipstream 
(references  6  through  8).  The  approach  used,  as  shown  in  Figure  10,  re¬ 
places  the  wing  by  a  system  of  discrete  small-span  horseshoe  vortices 
with  the  bound  vortex  elements  placed  along  the  1/ 4-chord  line  (assumed 
unswept).  The  wing  vortex  system  will  produce  velocity  perturbations 
which  in  general  violate  both  the  pressure  and  normal  velocity  boundary 
conditions  across  the  slipstream  interfaces.  To  correct  for  the  pressure 
boundary  condition,  a  system  of  horseshoe  vortex  elements  is  also  placed 
around  each  cylindrical  vortex  tube  representing  the  boundaries  of  the 
slipstream.  Violation  of  the  normal  velocity  boundary  condition  across 
the  slipstream  is  automatically  suppressed  by  using  a  reduced  velocity 
potential  ~  Kf  s  inside  the  slipstream,  as  done  by  Ribner  and  Ellis. 

The  tangential  flow  boundary  condition  is  satisfied  on  the  wing  at  the  3/4- 
chord  points  according  to  the  well-known  Weissinger  lifting  surface  theory, 
e.  g.  ,  Gray  and  Schenk  (reference  15).  However,  the  resultant  velocity 
includes  effects  of  slipstream  inclination  [Eqs.  (19)  and  (20)],  swirl 
( Eq.  (34)],  and  perturbations  from  the  slipstream  horseshoe  vortex  sys¬ 
tems,  in  addition  to  perturbations  from  the  wing  horseshoe  vortex  system 
as  in  reference  15.  The  pressure  boundary  condition  is  satisfied  across 
the  slipstream  in  a  linear  approximation,  as  originally  suggested  by 
Koning  (reference  1).  Satisfaction  of  the  wing  and  slipstream  boundary 
conditions  leads  to  a  system  of  linear  simultaneous  equations  whose  un¬ 
knowns  are  the  strengths  of  the  wing  and  slipstream  vortex  elements. 
Subsequent  solution  of  this  system  leads  to  the  span  loading  and  velocity 
field. 

A  further  simplification  has  been  to  take  all  trailing  vortex  elements  as 
straight  and  inclined  downward  to  the  free  stream  at  the  same  angle  6 
as  the  axis  of  the  vortex  tube.  It  is  recognized  that  a  somewhat  improved 
model  might  be  to  take  the  direction  of  the  trailing  vortices  as  coincident 
with  the  local  flow  angle  at  the  particular  spanwise  station.  Because  of 


the  additional  computational  complexities,  it  was  decided  to  forego  this 
remodeling  until  justified  by  comparison  with  experimental  data. 

In  the  following  sections,  expressions  are  derived  for  the  wing  and  slip¬ 
stream  boundary  conditions  in  terms  of  the  unknown  vortex  strengths. 
Results  of  the  computational  procedure  are  then  compared  with  test  data. 


BOUNDARY  CONDITIONS 


The  boundary  conditions  to  be  satisfied  are  a  constant  pressure  and  a  con¬ 
stant  normal  velocity  angle  across  the  slipstream,  and  tangential  flow  to 
the  effective  airfoil  mean  line  at  the  3/ 4-chord  position  of  the  wing. 

Slipstream  Boundary  Conditions 

The  pressure  boundary  condition,  from  the  incompressible  form  of  the 
Bernoulli  equation,  is  now 


“  !"|lls +  “sr’2* 


r 

L 


+  (w 

o 


+v> 


where  the  single  subscripted  velocity  components  are  for  zero  loading  on 
the  wing  [Eqs.  (19)  and  (20)],  whereas  quantities  with  a  second  subscript 
r  refer  to  perturbations  due  to  the  vortex  systems.  By  Eq.  (10),  and 
also  neglecting  second-order  perturbations  in  accordance  with  the  usual 
linearity  assumptions,  the  above  expression  becomes 


u  u_+v  v_  +  w  w 
s  sT  s  sT  s  sT 


=  u  u  _  +  v  v_  +  w  w  _ 

o  or  o  or  o  or 


This  condition  will  be  satisfied  on  the  slipstream  in  the  region  of  the  wing, 
for  which  the  first  terms  on  the  left-  and  right-hand  sides  may  be  as¬ 
sumed  dominant.  In  terms  of  the  perturbation  velocity  potentials  ^  8p 

and  ?0p  inside  and  outside  the  slipstream,  this  becomes 

us(d^sr/dx)-uo(d^or/dx)  (35) 


which  is  the  usual  linearized  pressure  boundary  condition  originally  intro¬ 
duced  by  Koning.*  Equation  (35)  should  be  satisfied  for  all  x  in  the  vicin¬ 
ity  of  the  wing,  and  at  all  azimuthal  positions  6  around  the  slipstream. 


♦  The  full  nonlinear  pressure  condition  has  been  satisfied  in  the  inclined 
actuator  disc  theory. 
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However,  following  Koning  and  others,  Eq.  (35)  is  integrated  from  x  =  -  ® 
to  x  =  co  giving 


us  ’’sr  **■  R>  9*  =  uo  ^or^®’  R*  9* 


(36) 


where  <p  _  (-  ®)  =  <p  (-  ®)  =  0.  This  allows  the  pressure  boundary  con- 

8 1  OX 

dition  to  be  satisfied  in  the  Trefftz  plane.  We  note,  however,  that  Eq.  (36) 
satisfies  the  linearized  boundary  condition  [  Eq.  (35)]  in  an  average  sense 
only,  and  is  therefore  less  exact. 


The  boundary  condition  on  the  normal  velocities  [Eq.  (1)]  may  be  readily 
generalized  to  include  the  velocity  perturbations  due  to  the  vortex  systems. 
Assuming  that  the  axial  perturbation  velocities  are  small  compared  to 
us  and  ,  and  utilizing  Eq.  (1),  we  obtain 


(37) 


which  can  be  satisfied  on  the  slipstream,  r  =  R,  for  all  x  and  8  by 
the  placement  of  a  suitable  distribution  of  sources  and  sinks  of 
strength  proportional  to  the  normal  velocity  jump.  Following  Ribner  and 
Ellis,  a  reduced  potential 


*'ar  =  (uo/us)  ”sr  (38) 

is  used  inside  the  slipstream,  so  that  the  normal  velocity  jump  is  sup¬ 
pressed.  In  terms  of  q> '  ,  the  pressure  boundary  condition  [  Eq.  (36)  ] 

now  becomes  sl 


*sr*00'  R'  9*  =  <u0/us)2  ^or*00’  R' 

In  terms  of  the  potential  jump  =  ?0p  -  ? '  p  =  r(6),  where  T(0)  is 

the  strength  of  the  horseshoe  vortex  element  on  the  slipstream  at  azimu¬ 
thal  angle  8,  the  pressure  condition  finally  becomes 

r(0)  =  (i  -p2)  <p oI-  (®,  R,  8)  (39) 

where  it  is  recalled  that  u  3  u  /  u  . 

o  s 

In  the  present  treatment,  Eq.  (39)  will  be  satisfied  at  each  of  M  slip¬ 
stream  control  points  of  azimuthal  angle  0m  as  shown  in  Figure  10. 

For  reasons  of  simplicity,  the  formulation  will  be  given  for  single-  and 
two-slipstream  configurations  only.  Generalization  to  configurations  with 
four  slipstreams  is  straightforward,  and  has  been  incorporated  into  the 
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machine  computational  procedure. 

Because  of  the  symmetry  conditions  discussed  previously,  Eq.  (39)  may 
be  satisfied  only  on  the  starboard  wing  panel,  and  0^  may  be  limited  to 

the  range  0  <  9^  <  it  (0  <  9^  <  tt/2  for  a  single- slipstream  configura¬ 
tion).  Thus,  we  are  led  to  a  series  of  M  equations  of  the  form 

i  2  f 

r<0  )  >  p  .  r.+  p  .  ..  r  +  p  r-p'r-p1  ,,r 

m  2  tt  mj  j  m(-j)  j  mn  n  mn  n  m(-n)  n 

m  =  1,2,...  M  (40) 

where  the  repeated  subscript  signifies  a  summation.  The  elements 

P  represent  the  contribution  to  the  potential  at  0  ,  just  outside  the 

mj  1  m  J 

slipstream,  arising  from  the  trailing  vortices  of  the  j’th  wing  horseshoe 

vortex  element  of  unit  strength  and  centered  about  y  =  y^  .  Pmj  *s  8^ven 

in  complex  notation  by 

P  .  =  arg  (Z.  +  h.  -  Z  )  -  arg  (Z.-h.-Z  )  (41) 

mj  6  J  J  m  J  J  m' 

Here  2  h.  is  the  vortex  spacing  on  the  wing,  Z.  =  Yj  »  and 
Z  =  y  +  R  exp  (i  0  )  ,  where  y  is  the  coordinate  of  the  slipstream 
axis.  The  elements  P^  jj  represent  the  contribution  to  the  potential  at 
0m  due  to  the  image  vortex  system  on  the  opposite  wing  panel,  and  are 
obtained  from  Eq.  (41)  through  the  relation 

P  .  ..  (Z  h.)  =  -  P  .  (-  Z  ,  -  h.)  (42) 

m(-j)  j,  j  mj  '  j  j 


The  elements  P  represent  the  contribution  to  the  potential  at  0 

mn  r  m 

from  a  pair  of  trailing  vortices  of  unit  strength  located  on  the  upper  half 

of  the  same  slipstream  at  azimuthal  angles  0n  t  A0n/ 2  ,  and  also  from 

their  images  on  the  lower  half  of  the  slipstream.  The  vortex  pair  at 
=  0^  is  included.  Because  of  the  symmetry  arising  when  the  wing 

passes  through  the  center  of  the  slipstream,  viz.,  TfG^)  =  -F(-0n),  it 

may  be  demonstrated  that  P  reduces  to 
’  mn 


where  6  is  the  Kronecker  delta  symbol  (see  Ribner,  reference  6). 
mn  ' 

The  elements  P^n  occur  only  for  multi- slipstream  configurations,  and 


19 


represent  contributi  ie  potential  from  the  image  vortex  system  on 

the  upper  half  of  the  earn  on  the  opposite  wing  panel.  ^3'mn  may  be 

written  as 


P' 

mn 


arg  R  exp  i  (it  -  0  +  A  0  /2)-R  exp  (i  0  )-  2y 

n  n '  m  ’  p 

arg  I”  R  exp  i  (it  -  0r  -  A  @n/ 2)  -  R  exp(i  0m>-  2  yp 


(44) 


Similarly,  the  P^  represent  contributions  from  the  lower  half  of  the 

slipstream  on  the  opposite  wing  panel,  and  may  be  obtained  from  Eq.  (44) 
through  the  relation 


P'  ,  ,  (0  . 

m(-n)  n 


A0n)  =  - 


P' 

mn 


(-  0 


n 


-  AV 


(45) 


Wing  Boundary  Condition 


The  condition  that  the  resultant  flow  velocity  is  tangent  to  the  effective 
wing  surface  at  the  i'th  control  point  along  the  3/ 4-chord  line  becomes, 
in  the  slipstream  coordinate  system, 


tan 


eff}i 


6 


(wR/u)., 


i  =  1,2,  .  I  (46) 


Here 

<Vi 


6  is  the  inclination  of  the  vortex  tube  representing  the  slipstream, 
is  the  resultant  Z  component  of  velocity,  and  (cr^)  .  is  the  effec¬ 


tive  angle  of  attack  of  the  wing  section  measured  with  respect  to  the  re¬ 
mote  free  stream  direction  ^  .  In  accordance  with  thin  airfoil  theory, 

we  define  as  the  angle  of  attack  of  an  equivalent  flat  plate  airfoil 

which  yields  the  same  lift  coefficient  as  a  cambered  and  flapped  section. 
Thus,  we  take 


a 


eff 


a 


-  °a  +  6f 


+  c  /2  TT 


(47) 


where  a  is  the  geometric  angle  of  attack  and  may  vary  along  the  span 

due  to  geometric  and  aeroelastic  twist,  «a  is  the  ideal  angle  of  attack 

of  the  section  (the  angle  at  which  the  pressure  at  the  leading  edge  of 
the  cambered  mean  line  is  finite),  is  the  section  flap  effectiveness 

parameter  at  zero  lift  coefficient,  6^  is  the  flap  deflection  angle,  and 

(3^  is  the  ideal  lift  coefficient,  i.e.,  the  lift  coefficient  at  a  =  a& .  The 

parameters  and  C  are  tabulated  in  reference  16  for  various  airfoil 
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mean  lines.  The  flap  effectiveness  is  also  given  in  reference  16  as 
a  function  of  the  flap/ wing  chord  ratio. 


For  control  points  located  outside  the  slipstreams,  we  take  at  the  i'th  sta¬ 
tion 

(u).  =  u 
1  o 


K 


'Vi  =  <w 


r>.*I 


(w  -  w  ), 

[  o  oo  k 


k=l, 2. 


+  w 


(48) 


where  (w^)  .  is  the  z  component  cf  the  velocity  induced  by  all  wing  and 

slipstream  horseshoe  vortex  systems  at  the  control  point,  and  wq  is 

the  z  velocity  component  due  to  the  k'th  inclined  slipstream  as  obtained 
from  Eq.  (19).  Thus,  the  multiple  slipstream  effects  on  the  basic  velocity 
field  are  approximated  by  superposition.  A  complete  theory  for  multiple 
inclined  actuator  discs  would  involve  additional  interference  velocities, 
inside  as  well  as  outside  the  slipstream,  and  has  as  yet  not  been  accom¬ 
plished. 


As  mentioned  in  the  previous  section,  a  reduced  perturbation  potential 
fp  =  (i(Pp  is  used  inside  a  slipstream.  In  order  that  the  form  of  Eq.  (46) 

remain  invariant  across  a  slipstream  boundary,  it  is  required  that  all 
other  velocity  components,  i.  e.  ,  in  addition  to  w^  ,  be  reduced  by  a  like 

factor.  Therefore,  for  control  points  inside  one  of  the  slipstreams,  say 
k  =  I  ,  we  take 


(u).  =  ll  u  =  u 
v  i  r  s  o 


(49) 
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ere  wsj  is  given  by  Eq.  (20)  and  Vg  by  Eq.  (34). 


All  quantities  in  Eqs.  (46)  through  (49)  can  be  regarded  as  known  at  a 
particular  control  station  except  for  (Wp,  ) .  Expressing  <Wp  )  .  in 

terms  of  the  T(yj )  and  F(9n)  gives 

<wr>i  =Sijr<Vj>4Sinr(°n>  (50' 

where  S.  and  S.  ^  are  influence  coefficients  which  give  the  z  velocity 
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component  induced  at  the  i'th  control  point  by  wing  and  slipstream  horse¬ 
shoe  vortex  elements  of  unit  strength  and  by  their  images.  Expressions 

for  S.  .  and  S.  nay  be  obtained  directly  from  the  Biot-S&vart  law, 

1  j  1  n  7  7 

e.  g.  ,  reference  1  7. 

VEIXDCITY  FIELD  AND  DOWN  WASH  ANGLES 


Once  the  strengths  of  the  vortex  elements  have  been  obtained,  the  result¬ 
ant  velocity  V  (x,  y,  z)  in  the  flow  field  may  be  obtained  in  a  straight- 
K 

torward  manner.  Outside  a  slipstream,  V  is  found  from 

K 


K 

=  V  (x, y,  z)  4  V 

k=  1 


(y.  z)  -  v 


4  v 


(51) 


where  Vj,  (x,  y,  z)  is  the  resultant  velocity  induced  at  the  point  (x,  y,  z)  in 

the  slipstream  coordinate  system  by  all  vortex  elements  on  the  wing  and 
slipstreams  and  is  obtained  from  the  Biot-Savart  law.  The  velocity  con¬ 
tributions  from  the  inclined  slipstream  are  obtained  by  superposition,  as 
indicated  by  the  summation  over  k. 

Inside  a  slipstream  the  induced  velocity  V^,  and  the  vortex  elements 
Tj  must  be  augmented  by  the  factor  ( 1  / p)  to  correspond  to  the  physical 
perturbation  velocity  and  vortex  strength.  We  therefore  obtain  for 
inside  the  f'th  slipstream 


K 

fR  =  <V.,=  ?  Vr(x.y,Z)^ 

k/f 


(52) 

Of  importance  in  the  estimation  of  the  static  stability  and  control  of 
V/STOL  as  well  as  of  conventional  aircraft  are  the  dynamic  pressure 
ratio  q/ q  ,  the  downwash  angle  t  ,  and  the  stability  parameter  dt/do. 

These  quantities  may  readily  be  found  from  V  For  convenience  in 

K 

presentation,  we  shall  express  these  parameters  in  free  stream  coordi¬ 
nates  (x1,  y'  ,  z'  )  rather  than  in  the  slipstream  coord  nates  (x,  y,  z) 
used  in  the  theoretical  development.  The  two  coordinate  systems  are  re¬ 
lated  by  a  simple  rotation  as  shown  in  Figure  1. 


vs  (y.  z) 


^v0(y,z) 
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The  dynamic  pressure  ratio  is  given  by 


q(x',y',z')/ q  =  V 


R 


Vr/V« 


(53) 


The  resultant  downwash  angle  €(x',y' ,  z1)  is  defined  as  the  angle  which 
the  projection  of  VR  in  the  plane  y' «  constant  makes  with  V^,  and  is 
given  by  an  expression  similar  to  Eq.  (33). 

The  parameter  d«/da  must  be  evaluated  in  a  body  fix?d  coordinate  sys¬ 
tem,  because  the  tail  is  assumed  to  rotate  rigidly  with  the  wing  as  the 
angle  of  attack  is  changed.  Thus,  dt/d  o  (x1  ,  y' ,  z')  is  evaluated  numeri¬ 
cally  from  the  expression 


di/da 


t(a,  x'  ,  y'  ,  z' )  -  t  (o  4  A  o,  x'  +  Ax' ,  y' 


A  z') 


/  A  a 


(54) 


where  Ax'  =  z'  A  O’  and  A  z'  -  -  x'  A  a  . 

RESULTANT  FORCE  AND  SPAN  LOADING 

Outside  the  slipstreams,  the  force  per  unit  span  n  may  be  found  from  the 
expression 


n  =  nQ  (y.)  =  p  T  (y.)  X  VR  q  (o,  y.,  o)  (55) 

where  x  =  z  =  0  corresponds  to  the  1/ 4-chord  line.  Inside  a  slipstream, 
because  of  the  augmentation  of  the  circulation  and  induced  velocities,  we 
take 


n  =  ns  (Yj)  =  P  (  r  (yj) /p)X  VRg  (0,  y.  ,  0)  (56) 

It  should  be  noted  that  n  is  the  resultant  of  the  local  lift  and  induced  drag 
per  unit  span,  because  induced  velocities  have  been  included  in  .  The 

direction  of  n  will  vary  across  the  span,  and  the  overall  lift  force  L 
and  drag  T5  may  be  obtained  by  integrating  the  components  of  n  perpen¬ 
dicular  to  and  parallel  to  V  across  the  span. 

The  span  loading  n(y)  must  be  continuous,  even  across  a  slipstream  in¬ 
terface,  because  of  pressure  continuity.  The  circulation  T(y^)  is  there¬ 
fore  discontinuous  across  the  slipstream  interface  (a  similar  discontinuity 
occurs  for  T(G  )  across  the  wing,  reference  7).  In  the  present  discrete 

vortex  formulation,  T(y.)  and  T(0  )  are  not  necessarily  con tinuous,  and 

jumps  in  their  values  will  automatically  occur  across  the  slipstream 
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interface  as  required  by  the  boundary  conditions.  This  is  in  contrast  with 
the  formulation  by  Ribner  and  Ellis  wherein  the  jumps  in  T  are  explicitly 
prescribed  and  can  unduly  complicate  the  analysis  for  multiple  slip¬ 
streams. 
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SAMPLE  CALCULATIONS  AND  COMPARISON  WITH  DATA 


Several  sample  calculations  have  been  carried  out  to  compare  the  present 
discrete  vortex  formulation  with  the  calculations  by  Ribner  and  Ellis  and 
with  experimental  data. 

A  comparison  of  the  calculated  span  load  distribution  with  that  from  refer¬ 
ence  8  is  given  in  Figure  11  for  a  wing  with  a  single  untilted  slipstream. 
According  to  Ribner  and  Ellis,  discontinuities  in  the  wing  vortex  strength 
are  required  across  the  slipstream  interface  r  =  R  such  that 

r  (R  +  o)  -  r  (R  -  o)  =  (i  -  p2 )  r  (R  +  o) 

whereas  a  similar  jump  in  slipstream  vortex  strength  is  required  across 
the  wing.  As  shown  in  Figure  11,  jumps  of  approximately  this  magnitude 
were  obtained  with  the  discrete  vortex  formulation. 

A  comparison  of  the  calculated  span  loadings  and  downwash  angle  with 
test  data  by  Stuper  (reference  4)  is  found  in  Figures  Id  through  16. 

Shown  in  Figures  Id  and  13  are  the  span  loading  and  downwash  angle 
distribution,  respectively,  for  a  wing  in  the  slipstream  of  an  untilted 
ducted  fan  with  straightener  vanes  to  remove  swirl  effects.  Good  agree¬ 
ment  is  found  except  at  a  =  ld*,  for  which  the  wing  has  probably  stalled. 

A  similar  comparison,  however,  for  a  wing  with  a  single  untilted  propel¬ 
ler,  is  given  in  Figures  14  and  15.  Because  the  present  theory  is  limited 
to  symmetrical  span  load  distributions,  slipstream  swirl  can  only  be  in¬ 
cluded  for  two-  and  four-propeller  configurations.  The  experimental  span 
loading  and  downwash  data  were  therefore  averaged  (right  and  left’  wing 
panels)  before  being  plotted  in  Figures  14  and  15  for  compa-ison  with  the¬ 
ory.  Once  again,  reasonable  agreement  is  found,  except  for  the  highest  a, 
where  stall  occurs,  and  except  for  some  differences  attributed  to  scatter. 

A  similar  span  loading  comparison,  for  a  single  propeller  inclined 
upward  6°  with  respect  to  the  wing  zero  lift  angle,  is  shown  in  Fig¬ 
ure  16.  Again,  average  experimental  data  were  used. 

Although  many  more  test  data  are  available  in  reference  4  for  compari¬ 
son  with  theory,  all  data  are  for  relatively  small  propellers  at  small 
inclination  angles.  Tests  of  comparable  detail  but  for  propeller  sizes 
and  tilt  angles  more  representative  of  V/STOL  configurations  and  opera¬ 
tional  conditions  are  required  to  substantiate  the  theory. 

This  requirement  has  been  partially  met  by  the  half-model  tests  of  refer¬ 
ence  14,  wherein  span  loading  data  were  obtained  by  testing  a  segmented 
wing  and  propeller  simulating  a  two-engine  V/STOL  aircraft  at  large  tilt 
angles.  A  comparison  between  the  present  theory,  including  effects  of 
slipstream  swirl,  and  representative  data  from  reference  14  is  given  in 
Figures  17  and  18.  Satisfactory  agreement  is  shown  at  angles  of  attack 
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and  propeller  inclination  angles  at  or  below  20°  for  which  wing  stall  ef¬ 
fects  and  wind  tunnel  wall  interference  effects  are  probably  small.  The 
test  data  at  larger  incidence  angles  and  for  thrust  coefficients  Tc' >  0.90 
have  not  been  included  in  the  comparison,  because  stall  and  wall  interfer¬ 
ence  effects  have  been  observed  in  the  data. 

Although  no  general  downwash  angle  survey  is  available  for  large  tilt 
angles,  i.  e.  ,  on  a  scale  comparable  to  that  of  reference  4  at  small  tilt 
angles,  it  is  still  possible  to  obtain  average  values  of  the  downwash  angle 
at  the  tail  location  f.  im  tail-on  and  tail-off  wind  tunnel  test  data.  Such 
data  (reference  18)  are  compared  with  the  theory  in  Figure  19.  The  test 
model  was  a  1/ ll-scale  version  of  the  four-engine  tilt-wing  configuration 
shown  in  Figures  36  through  50.  The  experimental  downwash  angle  was 
taken  as  the  tail  incidence  angle  with  respect  to  the  remote  free-stream 
direction  at  which  both  tail-on  and  tail-off  pitching  moment  coefficients 
were  equal,  and  represents  an  average  value  across  the  tail  span.  The 
theoretical  downwash  angle  was  calculated  at  the  location  of  the  horizontal 
tail  c/4-line  in  the  plane  y/  (b/2)  =  0.  325  (7 0°/i  of  the  tail  semispan).  Cor¬ 
rections  were  made  in  the  theory  for  the  vertical  di splacement  of  the  wing 
wake  and  of  the  slipstream  interface  by  shifting  these  boundaries  an 
amount  equal  to  the  calculated  displacement  of  the  wing  wake  (see  Design 
Cha  r  ts). 

As  shown  in  Figure  19,  satisfactory  agreement  was  found  for  c  1  10°. 

For  test  conditions  (a,  T  ",  &,)  which  resulted  in  downwash  angles  in  ex- 

O  C  I 

tess  of  10  ,  the  experimental  values  for  c  were  significantly  greater 
than  calculated  from  theory.  One  possible  explanation  for  this  disagree¬ 
ment  is  that  the  slipstream  and  wing  wake  are  deformed  in  such  a  manner 
that  their  downward  displacement  below  the  vertical  tail  is  appreciably 
reduced.  Thus,  it  is  known  that  a  completely  rolled-up  wake  is  displaced 
downward  by  only  2/  r  ^  of  the  displacement  of  a  plane  vortex  wake  from 
an  elliptically  loaded  wing  (reference  19).  Although  wake  roll-up  gener¬ 
ally  occurs  behind  the  tail  position  at  moderate  lift  coefficients  and  large 
aspect  ratios,  slipstream  effects  might  well  accelerate  the  roll-up  pro¬ 
cess,  because  of  the  low  effective  aspect  ratios  of  the  wing  segments  in¬ 
side  the  slipstream  (see  Introduction).  Similar  results  were  noted  by 
Heyson  and  Katzoff  (reference  20),  who  observed  a  very  rapid  roll-up  of 
the  rotor  slipstream.  They  found  that  the  rolled-up  vortices  had  been 
displaced  downward  only  about  one-half  as  far  as  the  calculated  momen¬ 
tum  wake. 

In  lieu  of  a  complete  wake  theory,  which  includes  slipstream  deformation 
and  wake  roll-up,  it  was  decided  to  correct  for  these  effects  by  using  only 
one-half  of  the  calculated  undeformed  wake  displacement  (or,  equivalent¬ 
ly,  to  shift  the  tail  downward  by  this  amount).  Substantially  improved 
agreement  with  the  data  from  reference  18  was  thereby  obtained,  as 
shown  in  Figure  20. 

Because  of  the  lack  of  sufficient  test  data  at  large  tilt  angles,  it  was  de¬ 
cided  to  forego  any  modification  of  the  theory  for  slipstream  deformation 
and  roll-up  effects  until  further  evaluation  and  comparison  with  test  data 
could  be  carried  out. 


DESIGN  CHARTS 


Knowledge  of  the  structural,  aerodynamic,  longitudinal  stability,  and  con¬ 
trol  characteristics  of  V/STOL  aircraft  requires  information  on  the  span 
load  distribution  and  on  the  wake  velocity  field.  In  order  to  supply  the 
aircraft  design  engineer  with  data  for  making  a  rational  estimate  of  these 
characteristics,  design  charts  have  been  prepared  for  the  span  loading 
C^'  (y),  downwash  angle  c  ,  stability  parameter  dc/do,  and  dynamic 

pressure  ratio  q/  q^  as  obtained  from  the  present  lifting  surface  theory 

(Figures  21  through  50). 

The  design  charts  were  prepared  for  twin-propeller  and  four-propeller 
tilt-wing  configurations  with  wing  and  flap  planforms  taken  similar  to  the 
Canadair  CL-84  and  the  Ling-Temco- Vought  XC-142a  aircraft,  respec¬ 
tively.  The  flight  conditions  and  aerodynamic  parameters  (e.  g.  ,  a,  a  , 

Tc’,  6^.,  J,  etc.)  were  chosen  as  representative  of  V/STOL  aircraft 

operations  in  the  hover,  transitional,  and  cruise  regimes.  Thus,  the  fig¬ 
ures  may  be  considered  an  extension  of  the  cruise  design  charts  by  Silver- 
stein  and  Katzoff  (reference  21)  to  conditions  wherein  inclined  slipstream 
effects  are  of  significance. 

The  data  for  the  charts  were  obtained  from  the  computer  program  des¬ 
cribed  in  Appendix  HI.  Running  time  was  approximately  30  seconds  per 
case.  The  data  for  c,  dt/da,  and  q/ q^  are  presented  in  the  form  of 

contour  plots  in  the  vertical  planes  y/(b/2)  =  0.40  and  0.325,  which  were 
assumed  representative  of  the  average  spanwise  horizontal  tail  stations 
for  two-  and  four-propeller  aircraft,  respectively.  Considerable  data 
smoothing  was  required,  because  the  discrete  vortex  formulation  gives 
spurious  results  at  points  near  the  singularities.  Although  all  computa¬ 
tions  are  carried  out  in  the  slipstream  coordinate  system  (x,  y,  z),  the 
contour  plots  are  presented  for  convenience  in  the  free  stream  coordinate 
system  (x1  ,y' ,z')  as  discussed  previously. 

It  is  noted  that  no  wake  roll-up  or  slipstream  distortion  effects  are  in¬ 
cluded.  In  this  respect,  the  charts  are  similar  to  those  in  reference  21. 
However,  because  of  the  large  angles  of  attack,  effects  of  the  vertical 
displacement  of  the  wake  have  been  approximated  by  shifting  the  contour 
plots  by  an  amount  Az(x,  y)  where 


x 


WR  (x,  y,  0)  d  x  . 


This  results  in  both  the  slipstream  boundaries  and  wing  vortex  wake  being 
displaced  downward  by  an  amount  equal  to  the  displacement  of  the  wing 
vortex  wake.  As  pointed  out  in  the  previous  section,  slipstream  deforma¬ 
tion  and  wake  roll-up  will  reduce  the  vertical  displacement  below  that 


27 


shown  in  the  design  charts.  In  lieu  of  a  calculation  procedure  which  con¬ 
siders  slipstream  deformation  and  roll-up,  the  present  charts  might  be 
corrected  for  these  effects  by  displacing  the  tail  position  by  an  amount 
equal  to  the  change  in  the  displacement  of  the  slipstream  or  wing  wake. 
The  physical  displacement  has  been  estimated  to  be  approximately  cne- 
half  the  calculated  displacement  Az  shown  on  the  charts  for  large 
tilt  angles. 

The  design  charts  include  the  effects  of  slipstream  rotation  with  viscous 
core  effects  reducing  the  swirl  velocities  to  zero  along  the  streamline 
passing  through  the  propeller  axis.  All  computations  were  made  with  the 
value  v  =  0.  10  ft^/  sec,  as  obtained  from  wind  tunnel  data,  with  no  scal¬ 
ing  for  aircraft  size.  The  effective  turbulent  origin  12  was  taken  as 
200  ft.  and  400  ft.  for  the  two-  and  four- s Upstream  configurations,  re¬ 
spectively. 

Jumps  in  the  quantities  «,  d  e/do,  and  q/q^_  occur  across  the  slip¬ 
stream  interface,  as  previously  noted  in  the  section  dealing  with  the  in¬ 
clined  actuator  disc  theory  (e.  g.  ,  Figure  5).  Values  of  d t/da  were  ob¬ 
tained  from  Eq.  (54)  using  an  incremental  A»  =  0.5°.  The  value  of 
d«/do  is  actually  undefined  at  the  slipstream  interface,  and  the  vaiues 
shown  are  limiting  values  obtained  by  extrapolation  from  inside  and  out¬ 
side  the  slipstream.  Viscous  effects  will  probably  play  a  role  in  smooth¬ 
ing  out  these  jumps.  Viscous  effects  are  also  known  to  modify  the  down- 
wash  angle  and  dynamic  pressure  contours  near  the  vortex  wake  of  the 
wing,  because  of  dynamic  pressure  losses  caused  by  skin  friction  on  the 
wing.  Corrections  for  the  latter  effect  are  given  in  reference  21,  and  re¬ 
quire  the  introduction  of  the  wing  profile  drag  coefficient  as  an  ad¬ 
ditional  parameter.  o 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A  theory  has  been  developed  for  ?n  inclined  actuator  disc  which  satisfies 
the  normal  velocity  and  nonlinear  pressure  boundary  conditions.  This 
theory  has  been  used  as  the  basis  for  the  development  of  a  wing  lifting 
surface  theory  with  inclined  slipstreams  for  application  to  V/STOL  air¬ 
craft  in  the  transitional  regime.  The  lifting  surface  theory  has  been  pro¬ 
grammed  on  the  CDC  3600  digital  computer.  Sample  calculations  have 
been  carried  out  to  compare  with  available  test  data  and  also  to  present 
design  charts  of  span  loading,  downwash  angle,  dt/da,  and  dynamic 
pressure  in  the  wake  for  typical  V/STOL  aircraft. 

Comparison  with  test  data  has  shown  that  the  theory  predicts  the  span 
loading  and  downwash  angle  reasonably  well  for  small  angles  of  attack  and 
small  propeller  tilt  angle.  However,  at  large  tilt  angles  the  theory  (with 
an  undeformed,  but  displaced,  slipstream  and  wake)  tends  to  predict  sig¬ 
nificantly  lower  downwash  angles  in  the  tail  region  than  observed  from  a 
single  set  of  test  data,  possibly  due  to  slipstream  deformation  and  wake 
roll-up  which  tend  to  reduce  the  downward  wake  displacement.  Use  of 
only  one-half  of  the  calculated  wake  displacement  gave  improved  agree¬ 
ment  with  measured  downwash  angles  at  these  conditions.  However,  the 
lack  of  sufficiently  detailed  test  data  for  conditions  representative  of 
V/STOL  transition  makes  evaluation  of  the  theory  and  assessment  of  its 
limitations  difficult. 

In  order  to  evaluate  the  theory,  a  test  program  is  required  for  a  wing- 
propeller  configuration  at  large  angies  of  attack  and  large  propeller  tilt 
angles  with  slotted  flaps  and  slats  to  prevent  separation.  Measurements 
should  include  span  load  distribution,  separate  overall  wing  and  propeller 
force  and  moment  coefficients,  and  a  survey  of  downwash  angle  and  dy¬ 
namic  pressure  inside  and  outside  the  slipstream.  Improvement  of  the 
theory  for  effects  of  slipstream  deformation,  wake  roll-up,  etc.,  could  be 
carried  out  following  experimental  justification. 
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Figure  2..  Control  Surface  for  Determination  of  Forces 
and  Power,  Inclined  Actuator  Disc  Theory. 
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Downwash  angle  e(  a>,  y,0)  ,  deg 


Figure  5.  Downwash  Angle  Across  Slipstream  Centerline, 
Inclined  Actuator  Disc  Theory. 
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Figure  8.  Comparison  W  ith  Slipstream  Rotation  Angle  Test 
Data  From  Reference  13. 
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Figure  9.  Comparison  With  Slipstream  Rotation  Angle  Test 
Pita  From  Reference  14. 
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Discrete  vortex  formulation 


Figure  11.  Comparison  of  Discrete  Vortex  Formulation  Loadings 
and  Vortex  Strength  With  Theory  of  Ribner  and  Ellis. 
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Figure  15.  Comparison  of  Downwash  Angle  With  Test  Data, 
Single  Untilted  Propeller  With  Wing. 
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Figure  17.  Comparison  of  Span  Load  Distribution  With  Test  Data: 

Two  Tilted  Propellers  With  Rotation,  T”  =  0.50. 
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Measured  downwash  an< 


Data  from  reference  18,  figs.  19  8  20. 

Measured  downwash  angle  is  tail  incidence  with 
respect  to  free  stream  direction  for  zero  tail 
moment  coefficient. 
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Two-Slipstream  Configuration  Design  Charts' 

T"  =0.5,  a  =  <?0°,  a  =  Z0°,  6.  =  0°,  J  0.  34. 
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Figure  25.  Two-Slipstream  Configuration  Design  Charts: 
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Two-Slipstream  Configuration  Design  Charts: 

T"  =  1.0,  a  -  80°,  a  =  90°,  6  =  0°,  J  =  0. 
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Two-Slipstream  Configuration  Design  Charts: 

T  "  =  0.  5,  a  =  40°.  o  =  40°,  6  ,  =  30°,  J  =  0.  34. 
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Figure  32.  Two-Slipstream  Configuration  Charts: 

T"  =  0.75,  a  =  40°,  a  =  40°,  6,  =  30°, 


J  =  0.  26. 
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Figure  33.  Two-Slipstream  Configuration  Charts: 
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Figure  35. 


Two-Slipstream  Configuration  Charts: 

T"  -  0.75,  o  -  50°,  a  .  60°,  b(  -  30°, 

c  P  1 


J  =  0.  26. 
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Four-Slipstream  Configuration  Design  Charts: 

T  "  =  0.75,  a  =  20°,  a  =20°,  b.  -  0°,  J  =  0.  37. 
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Figure  47.  Four-Slipstream  Configuration  Design  Charts: 

T"  =  0.75,  a  =  20°,  a  _  --  20°,  6r  =  55°,  J  =  0.26. 
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Figure  49.  Four-Slipstream  Configuration  Design  Charts: 

Tc"  =  0.  75,  a  =  40°,  a  =  40°,  6f  =  55°,  J  =  0.  26. 
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Figure  50.  Four-Slipstream  Configuration  Design  Charts: 

T  ”  =  0.75,  a  =  40°,  ap  =  30°,  6f  =  55°,  J  =  0.26. 
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APPENDIX  I 


EVALUATION  OF  INTEGRALS  ARISING  IN 
INCLINED  ACTUATOR  DISC  THEORY 


From  Eq.  (24)  we  have  defined 


I 

o 


/  9  z)  d  y  d  z 


Integrating  first  with  respect  to  z  , 
we  obtain 


and  recognizing  that  lim  <p  =  0, 

Z'**!  CD 


I 

o 


R 


“7  f  *0  |Zu'  dy 

IT  R 


where  z  is  the  value  of  z  ale  .g  the  upper  boundary  of  the  slipstream. 

u  2  2  2 
Making  use  of  Eq.  (18)  with  y  +  z  =  R  ,  we  obtain 


I 

o 


z  d  z 


where  a  =  2  ^/  ( 1  +  p) . 

No.  22a  of  reference  22), 


The  integration  may  be  evaluated  (e.  g.  ,  p.  86, 
and  we  obtain 


io  =  -0-n> 


Similarly,  the  integral  I  as  defined  in  Eq.  (24)  may  be  written 


I 
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Uu)  dy 


)  dy 


I 

o 


82 


where  we  recall  d<p  /  d  z  =  (p/w  )  (dq>  /3  z)  -  (1  -  p),  ar*8  we  ^ave  rnade 

s  00  s 

use  of  Eq.  (17)  evaluated  on  the  boundary. 

The  expression  for  Gq  may  be  defined  as 

Go  ■ -hfT 

TT  R  *'»' 

b2  o 

The  gradient  term  may  be  evaluated  from  Eq.  (18).  However,  it  is  more 
convenient  to  evaluate  y  <p  ^  directly  from  the  complex  potential  Fq. 
Thus,  in  polar  coordinates  we  obtain 

V^o2  =  (d  F()/dZ)  (dFj/dZ*)/wj 

B _ _ 

rS  Zr^  R  p.  cos  2  9 

Substituting  into  the  expression  for  Gq,  we  obtain 

— nr 

**  r  +  p  +  2  r  pcos  20 
1  o  r  r 

4  2 

The  integral  with  respect  to  0  is  2  -rr/  ( r  -  p  )  (reference  22).  Thus, 
G  becomes 


Evaluation  of  Gg  may  be  accomplished  in  a  similar  fashion  through  use 
of  the  complex  potential  for  Fg  [Equation  (15)]  and  through  use  of  the 

definitions  for  3<p  /  3z  and  <)<p  /  ">y.  This  gives 

G  =  G 
s  o 
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APPENDIX  II 


SLIPSTREAM  SWIRL  VELOCITY 


An  idealized  vortex  model  of  the  propeller  slipstream  (e.  g.  ,  reference  3) 
is  that  of  a  bound  vortex  of  constant  strength  T  rotating  with  the  blade. 
The  vortex  is  shed  from  the  blade  tip  and  moves  rearward,  because  of  the 
action  of  the  free  stream  and  because  of  self-induction  effects,  and  forms 
a  gpiral-like  trailing  vortex  system.  A  semi-infinite  vortex  of  equal 
strength  and  lying  along  the  slipstream  axis  must  be  inserted  into  the 
model  to  provide  continuity  with  the  bound  vortex.  The  axial  vortex  in¬ 
duces  swirl  velocities  inside  the  slipstream.  The  trailing  vortices  may 
be  shown  to  have  no  effect  on  swirl  inside  the  slipstream,  but  shield  the 
axial  vortex  outside  the  slipstream,  thereby  reducing  the  resultant  swirl 
velocities  to  zero  as  discussed  in  reference  3. 

Far  enough  behind  the  propeller,  the  trailing  vortex  spiral  may  be  as¬ 
sumed  infinite  in  extent.  The  axial  velocity  increment  induced  inside  the 
tube  [  see  also  Eq.  (3)]  is  then 

u  -  u  =  v  =  2N  r/(u  -t  u  ) 
so1  'os 

where  the  vortex  strength  per  unit  length  y  is  obtained  by  dividing  T  by 
the  number  of  spiral  loops  per  unit  length.  Here  N  is  the  number  of 
blade  revolutions  per  unit  time  and  (u  +  u  )  /2  is  the  average  axial 

05  2  2 

velocity  of  the  trailing  spi ral.  Since  AH/p  =  (1/2)  {u  -u  ),  we  obtain 


where  J  is  the  advance  ratio. 

Due  to  the  axial  vortex,  the  inviscid  swirl  velocity  Vq  will  be  taken  as 

v  q  /  VCD  =  (T/2irr  VJ  =  (J/2ir)  (Rp/r)  (AH/  qj 

which  is  the  first  factor  in  Eq.  (34).  The  second  factor  accounts  for  vis¬ 
cous  core  effects  (reference  12),  and  gives  a  zero  swirl  velocity  at  the 
center. 


84 


1 


APPENDIX  III 

FORTRAN  COMPUTER  PROGRAM 


GENERAL  DESCRIPTION 

An  extensive,  fairly  general  FORTRAN  computer  program  has  been  de¬ 
veloped  to  implement  the  methods  and  calculations  presented  in  the  body  of 
this  report.  While  written  for  and  run  only  on  the  CDC  3600  computer 
(32K  core),  adaptation  to  IBM  series  7040  and  7090  machines  has  been  an¬ 
ticipated  and  should  require  no  extensive  modifications  outside  of  plotting. 
(Plotting  could  not  easily  be  t ransformed .  ) 

The  general  organization  of  the  program  is  as  follows:  the  controlling 
program  calls  the  input  section  which  reads  and  stores  data  specifying  the 
problem  until  a  control  word  is  read  that  directs  a  specific  operation. 
Control  then  reverts  to  the  driving  program  which  supervises  the  desired 
calculation.  Most  numerical  results  are  output  as  they  are  derived; 
however,  with  the  plotting  option  active,  certain  results  are  retained  for 
further  processing.  If  the  plotting  option  is  on,  it  directs  the  smoothing, 
contouring,  and  plotting  of  these  retained  data.  The  main  program  then 
interrogates  the  input  routine  which  reads  and  stores  new  data  (if  any)  un¬ 
til  another  control  word  is  encountered.  This  processing  continues  until 
a  control  word  directs  the  program  to  stop. 

The  program  has  been  written  to  handle  any  wing-centered  configuration  of 
one,  two,  or  four  propellers,  provided  there  is  no  overlap  of  the  fully  con¬ 
tracted  slipstreams,  and  provided  the  span  loading  is  symmetrical  about 
the  airplane  centerline.  The  two-  and  four-propeller  cases  allow  for  slip¬ 
stream  rotation  as  in  Eq.  (34).  The  wing  has  of  necessity  a  straight  quar¬ 
ter-chord  line.  However,  the  chord  length,  twist,  flap-deflection  angle, 
flap  effectiveness,  ideal  angle  of  attack,  etc.,  may  vary  in  an  arbitrary 
manner  along  the  span.  These  spanwise  variables  are  input  in  a  flexible 
tabular  form  (to  be  described  in  detail),  and  linear  interpolation  is  used 
throughout.  The  plotting  option,  mentioned  earlier,  can  be  used  to  grapoi- 
cally  produce  wing  loading  curves  and  contour  plots  of  constant  downwash 
angle,  constant  dynamic  pressure,  and  constant  derivative  of  downwash 
angle  with  respect  to  angle  of  attack.  Considerable  hand  smoothing  of  the 
contour  plots  is  usually  required,  especially  near  slipstream  and  wake 
boundaries. 

To  complete  the  description  of  the  program,  the  input  format  and  options 
will  be  detailed.  A  sample  run  with  the  printed  input  and  output  has  been 
inc luded. 

INPUT  DESCRIPTION 

Input  to  the  digital  computer  program  is  read  first  as  an  alpha-numeric 
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card  and  checked  for  specific  control  words  in  columns  1  through  8.  If 
there  is  a  data  control  word  indicating  the  presence  of  a  data  list,  the 
card  is  re-read  using  the  format-free  NAMELIST  feature  of  FORTRAN  IV 
version  13.  {For  details  on  the  many  forms  data  values  can  take,  the  user 
is  referred  to  the  IBM  FORTRAN  IV  Reference  Manual.  )  If  the  card  has  a 
control  word  directing  a  specific  operation,  the  appropriate  action  is  tak¬ 
en.  If  there  is  no  dollar  sign  in  column  2,  the  card  image  (columns  3-80) 
is  taken  as  78  characters  of  BCD  title  information  to  serve  as  a  heading 
for  all  output  until  another  title  is  encountered. 

Th  e  program  is  organized  so  that  only  data  which  differ  from  case  to  case 
have  to  be  given.  Certain  default  values  are  provided  in  case  no  values  are 
given  explicitly.  For  example,  if  no  wing  twist  is  specified,  it  will  be  tak¬ 
en  as  zero  degrees.  The  default  values  for  the  various  parameters  are 
given  in  the  table  of  recognized  symbols,  normal  equivalents,  and  mean¬ 
ings  (Tables  I  -  III). 

The  problem  of  input  can  be  broken  into  three  parts: 

1.  Specification  of  the  wing  lifting  surface  and  propeller  configuration 

2.  Specification  of  the  remaining  quantities  (environment,  angle  schedule, 
and  numerical  control) 

3.  Direction  of  the  program  to  do  one  of  several  operations 
These  will  be  covered  in  the  following  sections. 

SPECIFICATION  OF  THE  AIRCRAFT  CONFIGURATION 


The  aircraft  is  treated  as  a  wing  with  a  straight  quarter-chord  line  and 
with  one,  two,  or  four  propellers  whose  fully  contracted  slipstreams  are 
centered  on  the  c/4  line.  Input  to  specify  this  lifting  surface  is  under  con¬ 
trol  of  a  data  control  word  A$WINGAA  (A  means  a  blank  character). 

The  span  of  the  wing  (B),  as  an  example,  could  be  given  as 

A$WINGA  A  B  =  1 00  $  END 

In  addition  to  the  span  (B),  certain  spanwise  functions  are  used  to  describe 
the  wing  geometrically  and  aerodynamics  lly.  The  chord  (CHORD),  the 
twist  (TWIST),  etc.,  are  all  considered  functions  of  y  (since  the  root  chord 
of  the  wing  is  in  the  xz  plane,  and  the  wing  extends  in  the  4y  and  -y  direc¬ 
tions).  The  spanwise  functions  are  stored  in  tabular  form,  and  linear  inter¬ 
polation  is  used  in  extracting  the  variation  from  these  tables.  Different 
sets  and  numbers  of  y  (Y)  values  may  be  associated  with  different  tables. 

A  single  entry  in  a  table  gives  a  constant  value  over  the  whole  span.  The 
maximum  number  of  entries  in  any  one  table  is  40,  and,  as  mentioned  ear¬ 
lier,  no  entry  will  result  in  a  default  value  in  some  cases  (see  Table  I).  It 
is  seen  that  the  input  form  is  arranged  to  conform  to  the  usual  section" 
description  of  a  wing,  but  this  is  not  required. 
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The  values  of  y  need  not  be  specified  with  each  table,  as  they  are  carried 
along  until  replaced.  Consequently,  the  entry  of  y  values  must  precede 
the  first  multiple  entry  spanwise  tabular  function.  The  y  values  given  on 
a  A$WINGAA  control  card  apply  to  the  following  A$WINGAA  data  cards 
until  a  new  set  of  Y's  is  read.  As  a  simple  example,  a  wing  with  a  60- 
foot  span,  a  5C  linear  washout,  a  10-foot  root  chord,  an  8-foot  tip  chord, 
and  a  constant  ideal  angle  of  attack  of  6.  67°  could  be  specified  as  follows: 

1  A$WINGAA  B  =  60,  ANGIDL(  1 )  =  6.  67  $END 

2.  A$WINGAA  Y  ( 1 )  =  0,  10,  r\\TST(l)=0,  -5  $F,ND 

3.  A$WING  AA  CHORD)  1 )  =  1 0,  8  $  END 

Several  points  are  illustrated  here.  All  numbers  and  symbols  in  the  data 
list  must  be  sepa-ated  by  either  "="  or  and  all  cards  are  terminated 

by  $END  (the  IBM  FORTRAN  IV  manual  describes  the  use  of  NAMELIST 
in  detail).  No  y  values  need  precede  the  spanwise  functions  with  only  a 
single  entry  but  must  be  given  before  the  double  entry  into  the  TWIST 
table.  A  table  may  be  continued  on  more  than  one  card  (columns  2-80) 
until  a  $END  is  read.  Since  the  y  value  table  is  destroyed  during  calcula¬ 
tion  of  a  case,  a  y  table  must  be  given  before  wing  data  for  any  case 
changing  these  values  (TWIST,  CHORD,  etc.).  It  is  seen  that  specifica¬ 
tion  of  wings  with  flaps,  taper,  etc.,  is  relatively  straightforward.  The 
reader  is  referred  to  Table  I  for  a  complete  list  of  wing  description  op¬ 
tions. 


SPECIFICATION  OF  THE  REMAINING  PARAMETERS  AND  DATA 

A  schedule  of  a  and  P  values  must  be  furnished  for  most  calculations. 
(The  program  uses  p^a^-a  to  specify  the  propeller  inclination.)  This  is 

accomplished  using  the  A$SCHEDA  control  card.  The  method  is  to  give 
the  starting  value,  the  final  value,  and  the  step  size  of  the  angle  in  de¬ 
grees.  For  example,  the  data  to  schedule  »  from  -5°  to  +10°  in  steps 
of  5°  and  P  from  -10°  to  +10°  in  steps  of  2°  would  read: 

A$SCHEDA  AZ=-5,  AMAX  =  10,  DA  =  5,  BZ  =  -10,  BMAX  =  1 0,  DB  =  2  $END 

The  symbols  recognized  following  a  A$SCHEDA  control  card  are  listed 
in  Table  III. 

All  remaining  data  are  input  under  the  direction  of  the  A$DATAAA  control 
word,  or  the  A$PLOTAA  word  used  to  control  the  plotted  output  options. 
All  the  symbols  recognized  following  a  A$DATAAA  word  are  listed  in 
Table  II,  and  all  those  recognized  following  a  A$PLOTAA  control  word  are 
listed  in  Table  III. 


87 


PROGRAM  CONTROL  CARDS 


There  are  several  control  cards  that  Instruct  the  program  to  perform  cer¬ 
tain  operations.  For  example,  A$STOPAA  will  terminate  the  program. 
The  input  section  of  the  program  reads  and  stores  data  until  a  control  card 
is  encountered  directing  some  operation  or  calculation.  When  that  action 
has  been  taken,  control  reverts  to  the  input  routine  which  continues  read¬ 
ing  cards  until  another  direction  is  indicated. 

The  A$RESETA  command  clears  all  of  the  data  from  the  storage  areas. 

It  is  useful  when  different  aircraft  are  to  be  processed  in  the  same  run. 
(The  results  of  A$RESETA  are  the  same  as  initial  startup  of  the  program 
with  exception  of  titles  and  case  number.  ) 

A$LOADAA  instructs  the  program  to  calculate  the  wing  loadings  only  for 
the  current  a  and  p  schedule.  Loadings  may  be  plotted  if  the  MPLOT 
flag  is  on. 

A$WAKEAA  instructs  the  program  to  calculate  loadings  and  to  calculate 
downwash  fields  and  their  derivatives  at  the  points  |  XMAP(i),  ZMAP(i)] 
for  each  YMAP(i)  A  whole  set  is  calculated  for  each  pair  (u,  (?)  in  the 
current  schedule.  Plotting  of  the  results  occurs  after  t-ach  case  if  se¬ 
lected  . 

A$STOPAA  terminates  the  run. 


SAMPLE  COMPUTER  INPUT  AND  OUTPUT 

Input  and  output  are  shown  on  pages  94-  96  for  the  four-slipstream  con¬ 
figuration  with  flaps  plotted  in  Figure  4  3. 
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TABLE  I.  RECOGNIZED  SYMBOLS  FOLLOWING  A 
A$WINGAA  CONTROL  WORD 

Program 

Symbol 

Normal 

Designation 

Default 

Value 

M  eaning 

Single  value  s 

ymbols 

B 

b 

*  * 

wing  span  in  feet 

[Independent  variable  for  spanwise  tabular 

functions 

*  Y  ( 1 ) 

y 

b/4 

the  current  values  of  y 
will  be  linked  to  any  table 
read  in 

Spanwise  tabular  functions 

♦  CMZERO(l) 

c 

m 

o 

0.  0 

wing  zero-lift  moment  co¬ 
efficient 

*  CHORDf  1 ) 

C 

*  * 

chord  length  in  feet 

*TWIST(1) 

°t 

0.  0 

twist  angle  in  degrees 

*CLIDL(1) 

c<a 

0.  0 

ideal  lift  coefficient 

*  FLPEFF(l) 

0  Or /  0  6j. 

1.0 

slope  of  the  curve  of  sec¬ 
tion  angle  of  attack  versus 
flap-deflection  angle  for 
zero  lift  (flap  effective¬ 
ness) 

*  FLPDEF(l) 

6f 

0.  0 

flap-deflection  angle  in 
degrees 

*  ANGIDL(  1 ) 

a 

a 

O 

O 

ideal  angle  of  attack 

Station  interface  specifications 

*  STA  TNS(  1 ) 

y/(b  /  2)  o, 

.05,  .1,  ... 

1.0  ratio  of  how  far  out  to 

place  horseshoe  interfaces 
relative  to  half  span 

*  These  variables  must  be  a  consecutive  list  of  from  1  to  40  values. 
*  *  Must  be  given. 


Program 

Normal 

Default 

Symbol 

Designation 

Value 

aning 

*  ARCS(l) 

0. 

l 

0,  30... 

180 

2  rops 

0,  15.  .  . 

90 

1  prop 

angles  in  degrees  at  which 
to  set  up  horseshoe  inter¬ 
faces  around  slipstream 


TABLE  II.  RECOGNIZED  SYMBOLS  FOLLOWING  A 
A$DATAAA  CONTROL  WORD 


Program 

Symbol 

Normal 

Designation 

De  fault 
Value 

Meaning 

VV  or  VELCTY  V 

CD 

must  be  given 
free  stream 

velocity  in  feet  per  second 

RHO 

P 

1.  0 

air  density  relative  to  a 
standard  day 

TCPP 

T  ” 

c 

0. 000001 

thrust  coefficient 

RADIUS 

R 

*  * 

radius  of  propeller  in  feet 

PROP 

VP 

0.  0 

spanwise  propeller  loca¬ 
tion  in  feet.  [  0.  0  means 
one  centrally  mounted  pro¬ 
peller,  nonzero  means  at 
least  two  propellers] 

PROP2 

'Vi 

0.  0 

0.  0  for  1-  and  2-propeller 
cases,  nonzero  gives  loca¬ 
tion  of  3rd  and  4th  propel¬ 
lers 

RPS 

N 

0.  0 

no.  of  revolutions  per  sec¬ 
ond  of  propeller  for  includ¬ 
ing  rotation  effects.  (Posi¬ 
tive  number  means  down¬ 
ward  rotation  inboard.) 

ETFVIS 

vt 

0.  1 

effective  viscosity 
(ft“/  sec) 

DPROP 

n 

100.  0 

effective  distance  ahead  of 
c/ 4  line  in  feet  to  locate 
propeller  for  purpose  of 
slipstream  rotation  effects 

STEP 

b/  10 

step  size  to  integrating 
trajectory  wake 

*  *  Must  be 

given . 
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TABLE  II  -  Continued 


Program  Normal  Default 

Symbol  Designation  Value  Meaning 


*XMAP(1) 

*  YMAP(l) 

*  ZMAP(l) 

PRINT 

various  digits  control 
certain  options  as  follows: 
unit  digits:  0,  1  track 

mapping 
y  stations 
only 

Z  track  all 
wing  and 
ring 
stations 
3  track  all 
wing,  but 
no  ring 
stations 
1-3  print 
wake 
profile 

tens  digit:  print  U  Y 
D(W  Y) 

hundreds  digit:  value  of 
MPLOT  kept 
here  internally 


x'  downwash  fields  are  cal¬ 

culated  at  each  lattice 

y'  point  given  by  these  co¬ 

ordinates.  Only  as  many 
z'  as  are  given  will  be  cal¬ 

culated. 

0  a  packed  integer  whose 


♦  These  variables  mav  be  given  in  a  list  of  up  to  10  variables. 


TABLE 

III.  RECOGNIZED  SYMBOLS  FOLLOWING  A$SCHEDA 
AND  a$PLOTAA  CONTROL  WORDS 

Program 

Normal 

Symbol 

De  signation 

Meaning 

A$SCHEDA  recognized  symbols: 

AZ 

Of 

mi  n 

control  scheduling  of  a  (wing 

AMAX 

a 

max 

angle  of  attack)  from  a  to 

6  min 

DA 

d  a 

a  in  steps  of  d  a  (degrees) 

max 

BZ 

Pmin 

’  control  scheduling  of  propeller  angle 

BMAX 

^max  ' 

with  respect  to  wing  from  p 

K  min 

DB 

dp 

to  p  ax  in  steps  of  d  p  (degrees) 

DEFAU LT 

SIZE 

must  be  given 
for  plotting 

total  height  of  contour  plot 

*WZC(1) 

no  downwash 
plot 

list  of  downwash  angles  to  contour 

*  DWZC(l) 

no  (It/  d  a  plot 

list  of  downwash  derivative  levels 
to  contour 

*  PC  ( 1 ) 

no  pressure  plot 

list  of  pressure  levels  to  contour 

M  PLOT 

no  plot 

mode  of  plotting 

0  no  plotting 

1-4  contours,  no  cross  plots 

5-8  both  contours  and  cross  plots 

9  cross  plots  only 
in  ranges  1-4  and  5-8  there  is 

1-4  point  interpolation  used  in 
producing  contours 

*  These  variables  may  be  given  in 

a  list  of  up  to  20  variables. 
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The  inclined  actuator  disc  analysis  is  combined  with  a  di  sc  r  ete  -  vo  rtex  Weissinger- 
type  lifting  surface  theory.  W  ing  -  prope  lie  r  combinations  at  arbitrary  wing  angle  of 
attack,  propeller  tilt  angle,  and  thrust  coefficient  are  considered.  Multiple  slip¬ 
stream  effects  including  slipstream  rotation  are  introduced. 

Agreement  between  theory  and  experiment  is  shown  to  be  satisfactory  for  small  slip¬ 
stream  inclination  angles.  However,  at  large  angles  the  theory  (with  an  undeformed, 
but  displaced,  slipstream  and  wake)  predicts  significantly  lower  downwash  angle 
the  tail  region  than  shown  by  the  test  data,  possibly  due  to  slipstream  deformation 
and  wake  roll-up.  Use  of  only  ine-half  the  calculated  wake  displacement  gave  im¬ 
proved  agreement  at  these  conditions.  However,  insufficient  data  are  available  for 
making  a  general  evaluation  of  the  theory  at  large  angles. 

Extensive  digital  computer  results  are  given  in  chart  form,  showing  span  loading, 
downwash  angle  <  ,  stability  parameter  d*/  d  u,  and  dynamic  pressure  at  aibitrary 
points  behind  the  wing  for  V/STOL  configurations  with  two  and  four  slipstreams. 
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